Enhanced inverse bremsstrahlung heating rates in a strong laser field 
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Test particle studies of electron scattering on ions, in an oscillatory electromagnetic field have 
shown that standard theoretical assumptions of small angle collisions and phase independent orbits 
are incorrect for electron trajectories with drift velocities smaller than quiver velocity amplitude. 
This leads to significant enhancement of the electron energy gain and the inverse bremsstrahlung 
heating rate in strong laser fields. Nonlinear processes such as Coulomb focusing and correlated 
collisions of electrons being brought back to the same ion by the oscillatory field are responsible for 
large angle, head-on scattering processes. The statistical importance of these trajectories has been 
examined for mono-energetic beam-like, Maxwellian and highly anisotropic electron distribution 
functions. A new scaling of the inverse bremsstrahlung heating rate with drift velocity and laser 
intensity is discussed. 

PACS numbers: 52.25.Dg,52.40.Nk,52.50. Jm 
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I. INTRODUCTION 

The description of electron-ion collisions in the pres- 
ence of a strong oscillatory electric field has been an 
essential part of every laser plasma interaction model. 
Processes such as collisional laser energy absorption, 
bremsstrahlung radiation emission and transport include 
electron scattering on ions in the presence of a high fre- 
quency electromagnetic field. Advances in laser tech- 
nology, particularly those related to the generation of 
ultra-short laser pulses and progress in inertial confine- 
ment fusion studies have challenged our understanding of 
scattering processes over a wide range of conditions and 
plasma parameters. Nonetheless, heating rate calcula- 
tions Iltillllllll[lllll[l3in strong laser fields have 
involved simplifications that are equivalent to a Born 
approximation for electron trajectories in the Coulomb 
fields. This approximation, although well-accepted in 
the physics community, should be revisited in view of 
recently discovered nonlinear processes such as harmonic 
generation and above threshold ionization [Til IT3. 
that have been explained as being due to strong modifi- 
cations of electron orbits in the laser field during electron- 
ion scattering events. Similar nonlinear modifications of 
particle trajectories have recently been discussed in nu- 
merical studies of inverse bremsstrahlung heating rates 
[13. Il5l|. Our paper continues this analysis; we examine 
classical electron trajectories that contribute to enhanced 
plasma heating and discuss their importance for cal- 
culations of effective collision frequencies with certain 
classes of electron distribution functions. We also sum- 
marize existing theories of inverse bremsstrahlung heat- 
ing and compare them with results of our test particle 
calculations. 

For the purpose of making comparisons, seminal works 
by Dawson and Oberman and Silin [3] and numerous 
subsequent publications on high field effects [1, S IE IE 



IE IE IE ^1 are classified into several broad categories: 
classical-mechanical analysis of binary electron-ion colli- 
sions in the dielectric approximation 0, 0, IE 13 ! kinetic 
theory which uses the Landau-type collision integral in 
the oscillating electron field (3]; ballistic model (or the 
model of instantaneous electron- ion collision) [3, l EJ l ^ ; 
" first principle" quantum-mechanical calculations [3, 1^ ; 
and, more recently, the quantum kinetic theory in the di- 
electric approximation ^5]. All theories consider the two 
body scattering problem as a starting point and their 
predictions lead to remarkably similar expressions which 
may differ only by a logarithmic factor. 

The agreement between these approaches follows from 
two common characteristics inherent in all theories, i.e. 
the Born approximation is used for the electron orbits 
in a Coulomb field and the cross-section is assumed to 
have a weak dependence on the laser field phase. The 
first approximation, with notable exception of the low 
frequency approximation |3|, has been employed in most 
quantum mechanical calc^ations where it applies to 
fast electrons with the velocity vq satisfying the condi- 
tion Ze^/hvo 1 01, where Ze is the ion charge. The 
Born approximation is equivalent to small angle scatter- 
ing and small momentum exchange in classical kinetic 
theory studies 0, |E • This corresponds to a straight 
line electron trajectory approximation, that is also a part 
of ideal plasma collision theories that lead to the Landau 
or Balescu-Guernsey-Lenard operators ^ simpli- 

fied description of electron orbits removes the sensitivity 
of scattering processes to the initial phase of an electric 
field. However, it has been demonstrated 0| that groups 
of electrons entering an ion interaction sphere at spe- 
cific phase can significantly modify the collisional cross- 
section and the heating rates. 

Irregular electron trajectories in the combined, high 
frequency field of a laser and a Coulomb field of an ion 
have been found by Wiesenfeld ^El- His numerical stud- 
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ies have identified stochastic trajectories for electrons 
with quiver velocity larger than the initial drift veloc- 
ity. This effect has been elaborate d up on in detail in 
recent studies by Fraiman et al. [T^ Il5l|. In addition to 
irregular trajectories that involve quasi-capture of elec- 
trons in complicated orbits encircling ions, a phase space 
analysis has revealed the existence of initial conditions 
corresponding to large angle correlated collisions with 
anomalously large energy transfer from the field to parti- 
cles. It was shown recently in Ref. (19(1 that the effect of 
correlated collisions exists not only in the Coulomb po- 
tential but also in short-range potentials in the presence 
of a strong laser field. 

There is a similarity between these processes and 
the interpretation of such atomic physics phenomena as 
multi-photon ionization and harmonic generation pro- 
posed by Corkum |^. It was suggested that free elec- 
trons that are created in the process of tunnel ionization 
can be brought back to the original atom by the laser field 
and undergo multiple correlated collisions. This process 
is sensitive to the phase of an electron in the oscilla- 
tory field and is responsible for processes such as high- 
harmonic generation, two-electron ejection, double ion- 
ization, etc. From the point of view of quantum mechan- 
ics, the high probability of secondary scattering was ex- 
plained by the phenomenon of Coulomb focusing jT^fl^ . 
This is a process by which the Coulomb attraction of the 
ion compensates for the natural dispersion of the electron 
wave packet. A classical analog involves the slowly mov- 
ing electron that is brought repeatedly into close proxim- 
ity of an ion as it performs large amplitude oscillations 
in the laser field. Each encounter results in a small de- 
flection of an electron trajectory toward the ion. These 
small changes accumulate into a deflection of an electron 
trajectory leading to a head-on collision producing large 
momentum change and causes an electron to leave the ion 
interaction sphere. An important question in the context 
of laser produced plasmas is determining under what cir- 
cumstances such trajectories are statistically significant 
and can alter macroscopic heating rates. 

Our paper examines effects of correlated collisions 
and irregular scattering trajectories on the inverse 
Bremsstrahlung heating rate. By using mechanisms such 
as the Coulomb focu sing [TTl ll^ . the parachute effect 
[T^ or quasi-capture |18| |. we analyze strongly modified 
electron orbits and classify them in terms of their ini- 
tial phase, impact parameters, and drift velocities. By 
means of test particle simulations we examine the sta- 
tistical importance of these trajectories. They indeed 
strongly modify the heating rates for all three different 
energy distribution functions used in the simulations, i.e. 
mono-energetic, Maxwellian and anisotropic electron dis- 
tribution functions. 

As we show in this paper, there are large discrepan- 
cies between existing analytical theories and test particle 
simulations of heating rates. They warrant further in- 
vestigations that may include molecular dynamics simu- 
lations IM, EH . Classical models of electron-ion interac- 



tions in molecular dynamics models include an effective 
potential to prevent the collapse of such system 

due to electrostatic attraction and properly account for 
quantum diffraction effects at short distances. We will 
also perform test particle simulations with this effective 
potential and will compare them with Coulomb potential 
scattering studies. 

The paper is organized as follows. In Sec. ^]we sum- 
marize existing theories of inverse bremsstrahlung heat- 
ing and present three representative expressions. Section 
mil presents a calculation model and gives examples of 
scattering orbits. Sectional describes the case of mono- 
energetic, beam-like electron distribution functions. We 
discuss different wave polarizations and field strengths. 
Heating rates for the case of a Maxwellian electron dis- 
tribution are calculated in Section Section IVII deals 
with anisotropic electron distributions that are encoun- 
tered in photo-ionized plasmas. Section IVIII contains a 
summary and conclusions. 



II. DESCRIPTION OF THE ENERGY 
EXCHANGE IN THE LASER FIELD 

In the introduction, we have identified several groups 
of analytical results describing electron-ion collisions in 
the presence of a homogeneous oscillating electric field. 
This division is a convenient characterization for various 
formalisms that are used to calculate the heating rate for 
electrons. All theories lead to remarkably similar heating 
rates, particularly for energetic electrons and small angle 
collisions, where electron scattering is well described by 
the Born approximation. 

The Dawson-Oberman model H H H H of the 
electron-ion correlation function is based on the classical 
Born approximation of binary collisions and on an ex- 
plicit average over random ion positions. The Coulomb 
field of the immobile ion of the charge Ze is considered 
as a weak perturbation to the electron motion which is 
comprised of a drift velocity vq along a straight trajec- 
tory and a quiver velocity v^; = eEi/muj in a laser field 
of amplitude E abd the frequency uj. Here — e and m are 
the charge and the mass of the electron. The Dawson- 
Oberman collision operator describes electron-ion scat- 
tering as an instantaneous event resulting in a small angle 
deflection of the electron trajectory and accounts for the 
screening effect of other electrons and the high frequency 
field in the dielectric approximation. The electron energy 
gain has been found in the weak field limit (jyj, and also 
for an arbitrary strength of the oscillatingfield Q with 
a Maxwellian distribution of electrons U13 or with an 
arbitrary distribution function Derivations by Pert 

7] and Shvets and Fisch arrive at similar expressions 
starting from the test particle, two body scattering prob- 
lem. The energy gain for a given electron in a plasma 



3 



with the ion density has the following form 0, : 



de 
di 



. fkvi 



k • vq 



OJ 



(1) 

where J/ is the Bessel function and the integration over 
the transferred momentum k involves the upper cut-off 
limit at k^ax = mv^j^Ze^ , which corresponds to the 
assumption of small angle scattering. 

Silin has calculated the high frequency nonlinear 
conductivity and the effective collision frequency for a 
fully ionized plasma by using a kinetic equation with a 
collision integral which accounts for small-angle scatter- 
ing. In this approach, the high frequency electric field 
defines the cut-off limit of the impact parameter and the 
effective collision velocity, but does not affect the colli- 
sion event itself. The Silin's collision frequency is in close 
agreement with the Dawson-Oberman model Q in both 
limits of weak and strong amplitude high frequency fields. 
A similar result for the electron energy gain has been ob- 
tained from the solution to the electron kinetic equation 
with the Landau collision operator after averaging the 
result with respect to the high frequency oscillations by 
Catto ^24!]. 

Several classical models of the effective collision fre- 
quency have been derived from explicit treatment of the 
dynamics of electron-ion scatterin g as suming instanta- 
neous and elastic interactions 0, IgL IIOI . We will identify 
them as the impact approximation Q 0| , or the ballis- 
tic model From the conservation of energy for such 
collisions, one finds that the change in the average elec- 
tron kinetic energy, which is related to the electron drift 
velocity, is proportional to the change in the electron mo- 
mentum. The energy gain is averaged with respect to the 
phase of laser field at the time of the electron- ion collision 
and expressed in terms of the transport cross section: 



dt 



A(v£ • Vq cos ^ + v\ cos^ 



Here, the Coulomb logarithm 
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depends on the full electron velocity v = vq -I- 
\ECO^<i> EE HI and prnax = V j uj and pjnm = 
max{2Ze^/mw^, Ti/mv} are the classical cut-offs at large 
and small impact parameters. According to Eq. ||2Jl, 
there is a singularity in the integrand if the quiver ve- 
locity, ve, approaches the drift velocity, vq. This singu- 
larity corresponds to orbits passing close to the ion and 
first appears for electrons launched in the direction of 
the electric field. Correlated collisions in presence of high 
frequency electric field have been observed in numerical 
simulations of particle trajectories [T^ll^ for the parallel 
launch, Vg || v^;. They remove this singularity and lead 
to an electron energy gain that is much larger than one 
that follows from Eq. ||2Jl, if < ^e- 



A quantum-mechanical description of electron-ion col- 
lisions in the Coulomb and laser fields employing the low 
frequency approximation has been developed by Bunkin 
and Fedorov ^ and KroU and Watson Agreement 
between the quantum mechanical and classical deriva- 
tions of the electron heating rate has been demonstrated 
by Ferrante et al. [2^. Within the theoretical frame- 
work of the Born approximation, this is a consequence of 
the well-known fact [23 that the electron scattering cross 
section in the Coulomb field is the same in classical and 
quantum-mechanical calculations. The energy gain av- 
eraged with respect to the Maxwellian electron velocity 
distribution function (2^] also shows a good agreement 
between classical and quantum results. The recent quan- 
tum theory of the inverse Bremsstrahlung heating rate 
PI, based on the dielectric approximation, compares well 
with classical results, provided the quantum cut-off pa- 
rameters are introduced into the classical kinetic theory. 

Quantum mechanical calculations _3, 28] describe the 
energy exchange rate between an electron and a laser field 
by using the cross sections derived by KroU and Watson 
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for electron-ion scattering. The result is given by 
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(4) 



where n and n' refer to the propagation directions of 
the electron before and after collision, respectively, ^ = 
2huj/mvQ, Vq is the initial drift velocity, and Imin = 1/C 
is the maximum number of emitted photons. Energy 
exchange in this equation is expressed in terms of the 
differences between the total absorption cross section (the 
terms with I > in the sum) and emission (the terms 
with / < 0) of Z photons. 



III. COMPUTATIONAL MODEL 

Our simulations of non-relativistic electron scattering 
in the field of an ion and in a uniform laser field are 
similar to studies of other authors 0, 0|. We solve 
Newton's equation for the motion of a test electron 



nir — — eEsinojt — WU{r) 



(5) 



where U{r) is the electrostatic potential describing the 
electron-ion interaction calculated along the electron tra- 
jectory, r{t). The electron velocity u = r is divided into 
two components: the quiver velocity, with amplitude v^;, 
and the drift velocity v and therefore u = v — v^; cos tot. 
In the absence of an oscillatory field and for the Coulomb 
potential, U{r) — Ze'^/r, the electron follows the Kepler 
orbit that is uniquely defined by its initial drift velocity 
Vq and the impact parameter p. The scattering process 
for the momentum exchange is described by the Ruther- 
ford cross-section, which corresponds to the surface area 
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of the interaction radius po = Ze^/mvQ multiplied by 
a logarithmic factor which accounts for small angle col- 
lisions. The limits in the logarithm are defined by the 
Debye screening length Xd = \/Te/4:ne'^ne at large dis- 
tances (here is the electron temperature) and pq. The 
small impact parameters correspond to large angle scat- 
tering events, including head-on electron-ion collisions, 
which contribute to the Rutherford cross section with 
finite probabilities. Similar properties describe electron- 
ion scattering in a weak oscillatory field ve/vq <C 1, the 
only exception being that the long distance cut-off is de- 
fined by the length vq/lu if it is smaller than the Debye 
length. 

The classical treatment of coUisional absorption and 
particle heating gives a result, which within logarithmic 
accuracy, agrees with quantum mechanical calculations 
in the Born approximation (5ll28j|. The exact agreement 
between the quantum scattering problem and the Kepler 
orbit solutions is well known and is a unique feature of 
the Coulomb potential Q^. In numerical studies based 
on classical mechanics such as molecular dynamics |23| 
or test particle simulations, the quantum diffraction ef- 
fect at short distances can be introduced by imposing a 
cut-off in the Coulomb potential at the electron de 
Broglie length, Xb = 2TTh/mv. This cut-off was first sug- 
gested by Uhlenbeck ^2^. It has allowed the successful 
application of classical theories to weakly degenerate sys- 
tems . For the simplest problem of electron scattering 
off a bare ion, no short range cut-off is necessary in view 
of the agreement between quantum and classical calcula- 
tions. No such general agreement exists, however, when 
the oscillatory field modifies the scattering orbits. For 
example, the quantum theory of inverse bremsstrahlung 
in the Born approximation 5] produces limiting expres- 
sions with logarithmic terms involving a cut-off at As. 
Also, there is no quantum theory going beyond the Born 
approximation '30l. Thus the problem of quantum cor- 
rections to the effective interaction potential is a difficult 
one and has no clear solution so far. 

For these reasons, the potential used in our simulations 
accounts for cut-offs at both short and large distances: 
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The two lengths, Xb and Xd, are considered as free pa- 
rameters that allow one to make contact with known lim- 
iting cases. Figure n compares the effective potential ^ 
with the Coulomb potential. The characteristic spatial 
scale of the problem is the distance rc = a/ Ze/E 
where the strength of the Coulomb field is the same of the 
oscillating field. Here, the Debye length was calculated 
for an electron temperature of I keV and an electron den- 
sity I0^° cm"-^. The de Broglie length, Xb — 2TTh,/mvE, 
was calculated at the quiver velocity ve corresponding 
to a laser intensity of 10^^ W/cm^ and a wavelength of 
0.25 /im. In our study, we examine, within the clas- 
sical formulation, scattering orbits that clearly display 
non-perturbative modifications to the Kepler trajecto- 



ries. We proceed with two dynamical models involving 
the Coulomb potential, as in , and the effective poten- 
tial jni with Xb calculated by using the quiver electron 
velocity. The upper cut-off at Xd is not important, be- 
cause the correlated collisions occur for relatively small 
impact parameters and do not demonstrate the logarith- 
mic divergence at large distances which is the character- 
istic feature of scattering in the Coulomb potential. 
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FIG. 1; Radial dependence of the effective potential U{r) @ 
normalized by Ze^ /rc in comparison with the bare Coulomb 
potential, rc/r. Examples of characteristic distances are 

shown in relation to rc = \/ Ze/E for Z = 10: the Debye 
screening length, Xd, was calculated for the electron temper- 
ature 1 keV and density 10^" cm~^; the electron de Broglie 
length. As, was calculated for the laser intensity 10^® W/cm^ 
and the wavelength 0.25 pm; the Bohr radius, as, and the 
classical electron radius, r^. 



Numerical solutions of the equations of motion, Eq. 
(O, require special attention regarding the choice of a 
suitable time step in computations of electron orbits near 
the scattering center. We have used the symplectic inte- 
gration algorithm |3lj |. which ensures a good numerical 
accuracy and an exact conservation of dynamical invari- 
ants. Even though the test particle problem described by 
Eq. (O appears simple, the construction of the statistical 
ensemble of results for a broad range of initial conditions 
for the average energy gain poses considerable challenges 
in terms of computer memory and time requirements. 
The overall geometry of our simulations is shown in Fig. 
121 Electrons were launched from the plane of incidence 
with a prescribed impact parameter p and an initial drift 
velocity vq which was always normal to the plane. The 
interaction sphere was defined by the screening length, 
Xd- The electron trajectories depend on the angle of vq 
with respect to the high frequency field polarization and 
on the field phase. 

The relevance of individual scattering trajectories in 
calculations of the average energy gain follows from the 
validity of the two-body collision approximation. This is 
also the basic assumption of all theoretical results used 
here for comparison. The two-body scattering model re- 
quires that the given electron can only interact with one 
particular ion during the scattering event. Therefore, the 



FIG. 2: Schematic electron orbits in the problem of electron- 
ion scattering in the laser field. Two typical electron tra- 
jectories are shown which originate at the launching plane. 
The sphere represents the long range cut-ofT of the effective 
potential. 



amplitude of electron oscillations, te = eE/muj^, in the 
laser field must be smaller than the average distance be- 
tween ions d = (3/47™^)^/'^. This condition, te < d, 
reads 



5.7 X lO^^-^A^ 



/on-/' < 1, 



(7) 



where the laser wavelength Ag is in /im, the intensity 
Iq is in W/cm^ and the ion density is in cm~^. For 
example, a laser pulse of intensity Iq ^ lO^'^ W/cm^ and 
wavelength Ao = 0.25 /zm satisfies the condition in a 
gaseous plasma of rii = 10^® cm~'^. 

The scattering trajectories that give rise to an enhance- 
ment of the electron energy gain |Q| are illustrated in 
Fig. O Figure compares three cases with the same 
initial conditions: a hyperbolic orbit in the absence of an 
osciUatory field (dotted hue), a trajectory in a Coulomb 
field and in an oscillatory field (dashed line) and a solu- 
tion to the equation of motion with an effective potential 
© and an oscillatory field (continuous line). In this ex- 
ample, the oscillatory field changes the trajectory and 
causes a head-on collision with the ion. This dramatic 
deflection is moderated by the quantum diffraction ef- 
fect accounted for in the effective potential ©. The sin- 
gle trajectory of Fig. corresponds to the case where 
the initial drift velocity is parallel to the direction of the 
oscillatory field. It illustrates the case of quasi-capture 
for the motion in the effective potential. The time de- 
pendent field reduces the instantaneous electron energy 
to a negative value allowing for a complicate bounded 
motion in the vicinity of an ion. The third example. 
Fig. 13;, displays two trajectories of electrons launched 
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FIG. 3: Examples of electron trajectories; (a) A hyperbolic 
orbit in the absence of the oscillatory field (dotted line) is 
compared with orbits for particles with identical initial con- 
ditions and parallel velocities with respect to the field polar- 
ization for the Coulomb potential (solid gray line) and for 
the effective potential (dashed line), (b) A complicated elec- 
tron orbit showing the quasi-trapping event - a long time 
trapping around the ion (solid line) for the parallel launch 
geometry, with the impact parameter equal to 2.5po and the 
initial drift velocity vq — O.Sve- This is compared with the 
hyperbolic orbit in the absence of the oscillatory field (dotted 
line), (c) Comparison of two similar trajectories of slightly 
different initial phases for the perpendicular launch geome- 
try. In this case the impact parameter is 15po and the initial 
drift velocity vo ~ O.Ive- The relevant Kepler orbit is shown 
as a dotted line. The remaining parameters for all examples 
are; Z = 10, the laser wavelength 0.25 pm and the intensity 
4.43- 10" W/cm^ 



in the transverse direction to the oscillatory field with a 
small difference in the initial phase. At first, both or- 
bits are very close to each other, yet after the collision 
they diverge dramatically. This sensitivity to phase and 
a strong modification of the scattering trajectory after an 
almost head-on collision, are manifestations of the effect 
of Coulomb focusing described in Refs. [Til IT^ Il3l| . 
Our calculations focus on the energy gain of a single 
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electron (e) for a given impact parameter averaged over 
the field phase and the rate of energy gain dt/dt for a 
single electron averaged over the impact parameters and 
the field phase. These quantities were calculated directly 
from the simulation results according to the following 
expressions: 

(^) = Y(^'o'--g>, j^^n,voJ{e)d'p, (8) 

where Vq and Vq are the electron drift velocities before 
and after the collision (cf. Fig. 12) and the angular brack- 
ets denote the average over the field phase. The electron 
energy gain from the oscillatory electric field © is com- 
pared with known analytical and numerical results. 

In the following three sections we describe the results 
of our test particle simulations of the electron energy 
gain for different initial electron distribution functions. 
We start with the beam-like mono-energetic electron dis- 
tribution function, followed by an isotropic Maxwellian 
distribution and finally with an anisotropic electron dis- 
tribution function. 



IV. ENERGY GAIN FOR A 
MONO-ENERGETIC ELECTRON BEAM 



Equation Q follows from Eq. after integration over 
the transfcrcd momentum k with the upper limit cut- 
off at kmax = mvo/h, which corresponds to quantum 
short distance cut-off at the de Broglie length calcu- 
lated at the initial electron drift velocity. The parameter 
Pmin = h/mvo- A morc detailed study of the quantum 
mechanical expression (|1(J|) can be found in Ref. [s^ . 
where it is shown that an enhancement of the total cross- 
section occurs for vq ^ ve and an oscillatory behavior is 
observed ior vq < ve- Note, that in the classical limit, 
<C 1, Eq. H1U|) gives an expression which is similar to 
Eq. ®. 

In the limit of a weak electric field, where ve/vq <C ^ <C 
1, one recovers from Eq. (|10|l the well known expression 

de ^ 47rn,Z^eS| / 2m«g \ 
dt mvQ \ huj J 

which accounts for the quantum cut-off at small distances 
and the high frequency cut-off, pmax — vo/^, at large 
distances. 

There is no easy way of simplifying theoretical results 
for the perpendicular polarization of the field. We use 
the original expressions ^ and Q for the comparisons 
with our numerical results. 



Here, we analyse the energy exchange between the os- 
cillatory field and electrons of a given initial velocity. 
The ensemble of different trajectories are constructed by 
varying the initial phase and the impact parameter. In 
particular, the parallel and perpendicular polarization of 
the field with respect to the initial direction of the drift 
velocity are considered. 



A. Analytical expressions 

The geometry of a parallel field polarization allows fur- 
ther simplification of the expressions for the electron en- 
ergy gain obtained by Dawson and Oberman, Eq. 
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and KroU and Watson, Eq. (0J: 
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B. Simulation results 

The results of numerical simulations are illustrated by 
Figs. 21|5l and 1^1 They display the average electron en- 
ergy gain, de/dt 0, as a function of initial drift velocity 
normalized to ve- Different vq/ve values along horizon- 
tal axis were achieved by changing the initial drift veloc- 
ity and keeping ve constant. In addition to vo/ve, the 
results also depend on the strength of the electric field. 
We characterize this dependence by using the parameter 
te/pe, where pE — Ze'^/mv'^ — pqVq/v'^ is the distance 
for a large angle scattering event for an electron with the 
velocity ve- It replaces po in the definition of scatter- 
ing cross-section for trajectories with vq/ve < 1- The 
parameter te/ Pe can be conveniently presented as 



teIpe = 3.4 X 10-20/3/2 a4/^ 



(12) 



(10) 



where the laser intensity is in W/cm^ and the wavelength 
is in jMn. It is large for typical conditions of present 
laser-plasma interaction experiments. However, the limit 
of small values of this parameter is also instructive to 
consider. 

Figure ^ corresponds to teIpe ^ 1 while Figs. EJ 
describe the high intensity limit, rE ^ Pe- As can be 
seen from Figs. 13 that are also drawn for te/pe ^ 1, 
scattering events of electrons experiencing large ampli- 
tude oscillations in the electric field are very sensitive 
to phase and are likely to result in large energy gains 
but only for a selected number of particles. Figure [51 
shows numerical results for the energy gain (jHJ with the 
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FIG. 4: The energy gain rate de/dt normalized to 
UiZ'^e'^ /mvE as a function of the normalized initial drift ve- 
locity for weak field values, rs ^ pE- The parameters are: 
Z = 10, the laser wavelength 0.25 pim and the intensity 
1.8 X 10^^ W/cm^. Dots describe numerical results for the 
parallel launch case (gray dots) and the perpendicular launch 
case (solid dots). They are compared with the KroU- Watson 
theory (dash-dot curves). Theoretical results showing nega- 
tive energy gains correspond to the parallel launch. 
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FIG. 6: The energy gain rate de/dt normalized to 
riiZ'^e*' /mvE for the parallel launch geometry. The param- 
eters are: Z = 10, the laser wavelength 0.25 /xm and the 
intensity is equal to 4.43 x lO^'^ W/cm^. Numerical results 
for the Coulomb potential are shown as large dots and for the 
effective potential as the solid curve; the Kroll- Watson ap- 
proximation - dash-dotted line, the Dawson-Oberman theory 
- dashed line, the classical approach - dotted line. 



increasing values of the energy gain for small initial ve- 
locities, ve which is in strong disagreement with all 
theoretical predictions. An agreement at large drift ve- 
locities follows from the applicability of the Born approx- 
imation to these trajectories. All curves in Figs. and 
El are well approximated by a classical result predicting 
the decrease of the energy gain asv^^. At small drift ve- 
locities, the calculated energy gain becomes larger. The 
increase is more dramatic for large fields, te/ Pe ^ 1, 
and for trajectories calculated with the Coulomb poten- 
tial (large dots) as compared to the effective potential ® 
(solid lines) in Figs. El and El 
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FIG. 5: The energy gain rate de/dt normalized to 
TiiZ^e'^/mvE as a function of the normalized initial drift veloc- 
ity for the perpendicular launch geometry. The parameters of 
the simulations are: Z = 10, the laser wavelength 0.25 /im and 
the oscillatory field intensity is equal to 4.43 x 10^^ W/cm^. 
Numerical results for the Coulomb potential shown as large 
dots and for the effective potential as the solid curve; the 
Kroll- Watson approximation - dash-dotted line, the Dawson- 
Oberman theory - dashed line, the classical approach - dotted 
line. 



Coulomb potential (large dots) and the effective poten- 
tial U (r) (EJ (continuous line) for the transverse launch 
(vq _L v^;) case. They are compared with results of the 
Dawson-Oberman theory ^ (dashed line), the classical 
model (dotted line), and the Kroll- Watson theory Q 
(dash-dotted line). The energy gain (jSJ for the parallel 
launch (vo|| v^;) is compared with results of three theo- 
retical models, (0), ©, and JTOl) in Fig. O 

Two main features characterize the numerical results 
in Figs. 21 El and El an agreement between theoretical 
and numerical results for high velocities, vq ^ ve, and 



C. Weak field limit 

We identify a weak field case by te/ PE ^ 1- Accord- 
ing to Eq. H12|l . this corresponds to relatively small in- 
tensities, /o < IO^Z^/^Aq**^^ W/cm^ where the laser 
wavelength is in /im. The energy gain rate is shown in 
Fig. 21 Electrons that are launched in either the parallel 
or perpendicular directions with respect to the oscilla- 
tory electric field lead to similar results. In both cases 
we observe an enhancement of the energy gain at small 
velocities. 

In order to interpret these numerical results, we first 
recall the electron scattering theory in the Coulomb field. 
The particle orbit is uniquely defined by the initial drift 
velocity, wq, and the impact parameter, p. The shortest 
distance along the Kepler trajecto ry that separates the 
electron from the ion is Tmin = vPo + P^ — Po- In es- 
timates which follow, we assume that the electron can 
gain a large amount of energy, on the order of mv'^/2, if 
it approaches the ion at least as close as p^. Then, by as- 
suming r^in = PEi we find the characteristic impact pa- 
rameter for these events pi„t = \/ Pe~^ '^PePo ~ \/'^PePo 
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for ve > vq. 

During the motion inside the interaction sphere, the 
electron oscillates many times with a relatively small am- 
plitude, te, and each time it undergoes a small angle de- 
flection from the original hyperbolic orbit. In an analogy 
with the Coulomb focusing process , each such deflec- 
tion brings the electron closer to the ion and eventually 
the electron undergoes a large angle scattering which ter- 
minates the interaction. As this trajectory modification 
involves many oscillation periods, it depends weakly on 
the initial phase. One could estimate the energy gain rate 
for the case ve > vq hy taking the product of the num- 
ber of scattering events per unit time niVQirpf^^ and the 
gained energy mv%. That gives the following estimate: 



<8> 



de/dt « 2'KniZ'^ / mvQ . 



(13) 



The 1/uo dependence results in large values of the en- 
ergy gain at small velocities in contrast to the previous 
theoretical results which do not depend on vq: de/dt ~ 
The correlation effects dominate 
the energy gain in the limit of small electron velocities, 
Vq < ve/A. 

The impact of this mechanism in the enhancement of 
the electron energy gain is greatly reduced for small in- 
tensities where ve/pe <C 1. Results of Fig. 0|have been 
obtained for the marginal case of r^; w pE- Still this gain 
enhancement is quite important as it is independent of 
the field polarization and the initial phase. That makes 
the weak field limit distinct from the regime te/pe ^ 1 
which is discussed in two next subsections. 



D. Strong field, perpendicular launch 

The interpretation of high intensity results {te/ Pe ^ 
1) is more complicated because of the strong dependence 
on the field phase. We first consider the energy gain of 
electrons launched transversely to the electric field di- 
rection. Fig. [SI Examples of relevant trajectories are 
presented in Fig. The energy gain is characterized 
by two impact parameters: px - in the direction the os- 
cillating field and py - in the direction of the initial drift 
velocity. As one can see in Fig. 13 the energy gain is 
an asymmetric function: it exists over a wide range of 
impact parameters px ^ te, while it is limited to much 
narrower range py ^ pint in the perpendicular direction. 

Modifications of electron orbits that lead to an en- 
hancement of energy gain also depends on the field phase 
(j). This is shown in Fig. |H1 which presents the depen- 
dence of the energy exchange on cf) for two initial drift 
velocities: the high velocity, vq = 3ve (solid line) and the 
slow initial motion, vq = O-Sve (dashed line). Assuming 
that incident electrons have a uniform phase distribu- 
tion, one concludes that only a small fraction of electrons 
~ Pe/{vo/uj) interact strongly with the ion while passing 
it within the distance of pE- A change in the electron 
energy can be negative or positive. This depends on the 
direction of electron motion, either towards or away from 
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FIG. 7: The average electron energy gain, (e) normalized to 
the Coulomb energy at the distance rc, Ze^ /rc = (Ze^ Ey^^ , 
as a function of impact parameters for the perpendicular 
launch geometry. The energy gain is averaged with respect to 
the field phases for vo = O.Sub. The parameters are: Z = 10, 
the laser wavelength 0.25 //m and the intensity 4.43 x 
W/cm^. The electric field points along the a;-axis. 



20 
15 
10 
5 



71/2 



71 



37x)2 



271 



-5 
-10 
-15 
-20 



FIG. 8: The electron energy gain, e normalized to (Z^E)^^^ 
as a function of the initial phase, for the impact parameters 
px = 2.5rc and py = 0. The perpendicular launch with vo — 
O.Sve (dashed lines) and 3ve (solid lines). Other parameters 
are: Z = 10, the laser wavelength 0.25 pm and the intensity 
4.43 X 10" W/cm^ 



ion when it enters the sphere of radius pe- For fast elec- 
trons, the energy loss for particular phases turns out to 
be exactly the same as the gain for other phases with no 
net change to the Born approximation results. For slow 
particles (dashed line in Fig. IHJ the reduction of electron 
energy can result in quasi-trapped trajectories. This mo- 
tion is irregular and produces a net energy gain when 
the electron finally escapes the ion. This effect is more 
important for slower electrons as it leads to increasing 
energy gain at small drift velocities as shown in Fig. [S] 

An estimate of the energy gain rate de/dt involves the 
characteristic cross section, a = irpxPy^ and the char- 
acteristic energy gain mv'^. The impact parameter, px, 



9 



along the field direction can be found by assuming that 
the upper bound for the closest approach and large angle 
scattering is te, that is, te — \/ Pq + Px ^ Po leading to 
the estimate px — rE\/l + 2po/rE- As discussed above, 
one have also to account for the fact, that only a small 
fraction of electrons, ~ pE/{vint/'-^), undergo collisions 
with this impact parameter. Here, the effective drift ve- 
locity Vint = vo vl + 2po/rB is calculated at the distance 
rE from the ion. By including these results in the expres- 
sion for the energy gain, one obtains 



de 
dt 



T:n,vo{2poPEy^^rEil + 2po/rB)^/^ x 

r72 4 

Pbcj TTUiZ e 

vq{1 + 2po/rEy^^ mvQ 



(14) 



This estimate is very similar to Eq. for the energy 
gain in the weak field limit. The I/vq dependence ex- 
plains the growth of the energy gain at small velocities. 
These are the large angle collisions of selected electrons 
that contribute to this enhancement. This is moderated 
by quantum diffraction effects, which are accounted for 
in the expression for the effective potential (solid line in 
Fig. [SJ if the interaction radius pint becomes smaller 
than the de Broglie length. 



E. Strong field, parallel launch 

The perpendicular launch considered above represents 
a general case - electrons entering the interaction sphere 
at arbitrary angles with respect to the polarization vector 
of the oscillatory field are gaining on average the same 
amount of energy as those launched in the perpendicu- 
lar direction. Only the electrons propagating parallel to 
the polarization vector make an exception. As shown in 
Fig. these electrons are gaining much larger energies at 
small velocities as compared to the previous cases. The 
explanation of this anomalous energy gain is related to 
Coulomb focusing JJ^ or the " parachute effect" 14, . 
Electrons launched in the parallel direction can return to 
the same ion several times. Each time, as the electron 
passes the ion, its trajectory is defiected towards the ion. 
The effective impact parameter decreases until a head- 
on collision terminates the interaction. The number of 
possible returns scales as ve/vq. 

This scattering geometry has been analyzed and de- 
scribed in detail by Fraiman et al. 0, 0| . A typical de- 
pendence of the energy gain averaged over the field phase 
(e) © on the impact parameter for the parallel launch is 
shown in Fig. [Sjfor vq = O.Sve- It displays three maxima 
in the energy exchange corresponding to three electron 
returns. The second curve (solid line) in Fig. I^jhas been 
obtained with the effective potential 10. The amount of 
energy gained by electrons depends on the short range 
behavior of the interaction potential at r ^ po and is re- 
duced when close distance interactions are cut-off (in the 
present example Xb > Po)- Figure [TUI illustrates the sen- 
sitivity of the electron energy gain to the field phase and 



the relative importance of head-on collisions. The latter 
are represented by the large maximum near = 37r/2. 
The effect of quasi-capture and random trajectories do 
not strongly contribute to the energy gain in this case. 




FIG. 9: The average electron energy gain, (e), normalized 
as a function of the impact parameter for the 
parallel launch. Numerical results are averaged with respect 
to laser field phase for vo — 0.3ve (dashed line). Other pa- 
rameters are: Z = 10, the laser wavelength 0.25 pm and the 
intensity 4.43 x 10^^ W/cm^. Solid lines are obtained from 
simulations with the effective potential. 
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FIG. 10: Dependence of the electron energy gain, e, normal- 
ized to {Ze^Ey^'^ on the initial phase for the impact param- 
eter 2.5rc for the parallel launch and for vo = 0.3«e. Other 
parameters are: Z = 10, the laser wavelength 0.25 ^m and 
the intensity 4.43 x 10^® W/cm^. Dotted lines correspond to 
simulations with the effective potential. 

An estimate of the heating rate can be made by fol- 
lowing the arguments of Fraiman et al. 14, I^j . When 
passing by an ion, an electron with an impact parame- 
ter p is scattered by a small angle S9 « 2pe/p in the 
Coulomb potential. This defiection changes the origi- 
nal orbit and the electron acquires a transverse velocity 
component v± = 69 ve- That electron will arrive into 
the head-on collision with the ion after half a field pe- 
riod if p = TTVj_/u}. From the above expression one finds 
p = \/2pETE = ^J2rc- The electrons with larger impact 
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parameters will also be attracted to the ion but only af- 
ter several oscillations. Since the number of possible re- 
turns is ve/vq, the characteristic impact parameter for 
this strong interaction is estimated as tcve/vo. Then 
the average electron energy gain for the parallel launch 
can be evaluated as 



dt 



PE 



— « 2TmiVo—r-'(j^mv% 
at rc 



2'Kni I — 



Pe 



1/2 



(15) 

where we have also accounted for the fact that only 
fraction of electrons proportional to pe/tc will undergo 
strong interactions (cf. Ref. Although the formula 

(|15|l demonstrates the same dependence on the electron 
drift velocity as ((1^ . the energy gain is (teI PeY^"^ ~> 1 
times bigger due to the parachute effect discussed above. 
This Coulomb attraction due to multiple collisions is ef- 
fective, however, only for sm all angl es between Vq and 
the electric field vector, < ^ Pe/te, where the velocity 
v±^ has a well defined sign. For larger angles, the electron- 
ion interaction reduces to a single collision event, as pre- 
sented previously for the case of perpendicular launch. 



limit, ve VTe, the effective collision frequency coin- 
cides with the electron-ion transport collision frequency 
Vei — 4-\/27rZ^e^niA/3m^'y|,g, where the Coulomb loga- 
rithm A is defined by Eq. Q. In the strong field limit, 
ve ^ VTe, i^eff dccrcascs as Vj^^ and there is a simple 



expression from Ref. |33l l34| which provides a very good 
interpolation for all existing analytical theories: 



(18) 
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V. HEATING RATE FOR AN ISOTROPIC 
ELECTRON DISTRIBUTION 

The energy gain rates H14|) and (|15|l calculated for 
mono-energetic electrons can be used for evaluating heat- 
ing rates for various electron distribution functions. Here 
we consider the common case of an isotropic Maxwellian 
electron distribution function which is established due to 
electron-electron and electron-ion collisions. Before ana- 
lyzing the statistical significance of correlated collisions, 
we recall the basic elements of the classical theories. 

By taking an average with respect to a Maxwellian dis- 
tribution function with the temperature Tg, the Dawson- 
Oberman theory Q gives the following formula for the 
heating rate per electron: 



dTf. , — 



4, ,2 +°° 



mw 



Te 



1=1 



-^Ii{q)exp 



(16) 

where Ii{x) — Jf{y)dy, vtc — \/TeJrn is the electron 
thermal velocity, and the parameter — 2Tujj /Tg defines 
the cut-off at small distances. The quantum mechanical 
theory gives the following heating rate 



dt 



flUJVTe 



I sinh [ 



—h{q) X 



exp 



' 2q^T, 



q^h^LU^ 



(17) 



One can define the effective collision frequency by the re- 
lation: dTe/dt = iyeffmv'^/2 and use Eqs. ifTCjl and lfT7|) 
to evaluate it. It is well-known that in the weak field 



dT/dt 
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FIG. 11: Dependence of the effective collision frequency 
normalized to the electron-ion collision frequency — 
4\/27rniZ^e''/3m^ii|ng (a) and the absorption rate per elec- 
, tron dTe/dt normalized to VeiTe (b) on ve/vtb for a plasma 
with Z = 10 and for the Maxwellian electron distribution 
function with Te = 200 eV (black dots) and 500 eV (grey 
dots). The continuous lines following large dots are calcu- 
lated from the formula II19II . Numerical results are compared 
with the theoretical results for Te = 200 eV: the KroU- Watson 
approximation - dash-dotted line, the Dawson-Oberman ap- 
proximation - dashed line, the classical approach - dotted 
line, and the expression l|18|l - the fine solid line. The laser 
wavelength is 0.25 fim. 

Figure [TTk. compares the variation of effective collision 
frequency with ve/vtc found from our simulations and 
from theoretical formulas presented above. The numeri- 
cal heating rates were calculated by the direct average of 
the energy gain rate over the isotropic Maxwellian distri- 
bution function. Numerical results have been obtained 
for the cases when the electron temperature is 200 eV 
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and also for 500 eV. All models and numerical simula- 
tions generally agree for ^ f Te • Some discrepancies at 
small values oi ve follow from differences between theo- 
retical approximations, in particular, from the treatment 
of short range interactions. 

The most dramatic difference between the known theo- 
ries and our results is illustrated in Fig. Illb which shows 
that the calculated heating rate reaches a constant val- 
ues in the large field hmit, ve > vtb- Our results for the 
effective collision frequency can be approximated by the 
following formula: 



.(1 



0.3 4/4, 



(19) 



This expression gives a 1 /vtc dependence in the heating 
rate for the limit of large values of ve/vte- This agrees 
with the conclusions of Ref. Such behavior agrees 

with the 1 /vo dependence of the energy gain rate for the 
perpendicular electron launch discussed in sections llV (]l 
and IIVDI The contribution of electrons with parallel 
velocities is not evident in Eq. (|19|l because small angles, 
^ ^ \J PeI^e^ make a small statistical contribution to the 
average quantity for the isotropic distribution function. 

We recall that the enhancement of the heating rate 
due to the correlation effect takes place in a sufficiently 
rarified plasma where the electron quiver radius is 
smaller than the inter-particle distance d. The condition 
Q ensures this limitation. 



VI. 



HEATING RATES FOR AN ANISOTROPIC 
ELECTRON DISTRIBUTION 



It was shown in Sec. IIV El that slow electrons (wq <'Ve) 
propagating along the oscillatory field are the main con- 
tributors to the enhanced heating rates. However this 
effect is confined to small angles and it disappears if 
electrons are distributed isotropically. Therefore for typ- 
ical laser plasma interaction conditions, the heating rate 
might depend dramatically on the number of slow elec- 
trons and on the anisotropy of actual electron distribu- 
tion function. In particular, this effect could be impor- 
tant for plasmas produced by the ionization of gase s b y 
high intensity laser pulses. According to Refs. |35L ISq. 
the velocity distribution function of the photoelectrons 
displays a strong anisotropy along the direction of the 
linearly polarized laser beam, if tunnel ionization is the 
dominant mechanism. The following simple analytical 
expression captures the main features of electrons pro- 
duced by the tunnel ionization process |3^ : 

Fo(v) = ne(m/2^r,)i/25(vi) cxp(-mt;f /2re) (20) 

where v\\ is the electron velocity in the direction of laser 
polarization and the effective electron temperature Te 
may be a function of the ionization potential and the 
laser field amplitude [3^ . In our analysis we treat Tg as 
a free parameter. 



We compare results of various analytical models with 
our simulations after averaging them with respect to the 
distribution function (|20|) . This averaging reduces to a 
simple integration of expressions for the electron energy 
gain in the parallel geometry lO and H10() with a one- 
dimensional Maxwellian velocity distribution. The main 
contribution to this integral comes from slow electrons. 
In the case where ve < VTe, all classical theories give 
similar results that depend on the cut-off parameters. 
For the quantum approach, the heating rate per electron 
derived from Eq. I|10|) has the following form: 



dt 



+ 00 

l^^E \ rfuexp 



^E „,2 



2T, 



u X 



dy 



{yu + ^Eim^' 



(21) 



n — y/u^+^jsl 

where ^e = 2?iLu/mv'^. 
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FIG. 12: Dependence of the absorption rate per electron 
dTe/dt normalized to feiTe on the electron quiver velocity 
for a laser wavelength 0.25 /im and for an anisotropic dis- 
tribution function (I20I I with Te = 200 eV: numerical results 
- dots, the KroU-Watson approximation - dash-dotted line, 
the Dawson-Oberman approach - dashed line, the classical 
approach - dotted line, expression (1221 - solid line. 

Figure IT^ shows the intensity dependence of the heat- 
ing rates that are obtained from our numerical model and 
from various analytical expressions. The striking discrep- 
ancy between the analytical results and our numerical 
simulations is a consequence of the large number of slow 
electrons in the electron distribution function (|20|l prop- 
agating parallel to the polarization direction. Contrary 
to the case of a three dimensional Maxwellian distribu- 
tion function, expression H20I) contains a finite number 
of electrons at zero velocity and therefore the heating 
rate diverges as a logarithm of Vmin , if one averages the 
energy gain (|15|l with the one-dimensional electron dis- 
tribution function. This divergence can be resolved by 
accounting for the fact that the maximum impact pa- 
rameter pmax ~ re Ve/ vq cannot be greater than the 
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inter-particle distance d. This leads to the following es- 
timate Vjnin ~ VETc/d. Then one finds the following 
expression for the heating rate per electron: 



dt mvTe 



PE 



1/2 



In d^ 



This expression is proportional to v^J.'^ . It approximates 
well numerical results as shown with the solid line in Fig. 

[El 

Photo-ionized plasmas with anisotropic electron dis- 
tributions provide an example of a case where enhanced 
heating rates are obtained due to correlated collisions 
as described in our paper. It is clear that the partic- 
ular form H20() of the electron distribution function can 
be quickly altered due to rapid heating. However, in- 
verse bremsstrahlung heating can potentially support a 
certain level of anisotropy in the electron distribution 
function. More work needs to be done in developing the 
self-consistent electron distribution for this case and in 
comparing this model with experiments. 



T, 



Ze^TE 



(22) 



VII. CONCLUSIONS 

We have performed numerical simulations based on the 
test particle model which involves classical electron tra- 
jectories in the electric field of an ion and a homogeneous, 
high frequency laser field. We have confirmed recent re- 
sults by Fraiman et al. regarding the electron orbits 
that closely approach an ion several times during each 
scattering event and subsequently lead to a large angle 
deflection and to a large enhancement in the electron en- 
ergy. These correlated collisions have been ignored in 
preceding theories which significantly underestimate the 
electron heating rates by neglecting large angle scatter- 
ing and the important field phase dependence of electron 
trajectories. The inverse-bremsstrahlung heating rates 
have been evaluated for different electron distribution 
functions, including: mono-energetic beam-like electrons, 
Maxwellian and highly anisotropic distribution relevant 
to photo-ionized gases. The energy gain for electrons 
with small drift velocities and for large amplitude elec- 
tric fields, vo < ve, increases as I/vq and exceeds the 
previously known results if vq < ve/^- This behavior 
contradicts results of all analytical theories and it is espe- 
cially important for electrons with initial velocities that 



are parallel to the field amplitude. When the quantum 
short-distance cut-off is introduced through an effective 
potential, the energy gains demonstrate a similar depen- 
dence on the particle drift velocity but a smaller mag- 
nitude of enhancement is obtained if the cut-off length, 
Xb — 2'Kh/mvE, is larger than the characteristic radius 
of the Coulomb interaction pc = Ze^/Tg. 

For isotropic Maxwellian distribution functions, the 
heating rate in the limit of > vtc is independent 
on the laser intensity. It can be approximated by the 
interpolation expression (|19|l for the effective collision 
frequency. For highly anisotropic electron distribution 
functions such as those produced by photo-ionization, the 
new effect of large angle correlated scattering dominates 

3 /2 

the heating rate at high laser intensities. It leads to v-^ 
dependence of the heating rate. 

Our numerical results involve calculations of the elec- 
tron scattering on a single ion. They are relevant to re- 
alistic plasma conditions provided the two-body collision 
approximation r^; < o? Q is satisfied. This imposes an 
upper limit on the laser intensity for a given ion density 
and typically holds for very underdense plasmas. Also, 
we have implicitly assumed that the plasma is weakly 
coupled and that the average distance between ions, d is 
smaller than the screening length, Xe- 

The enhancement of the inverse bremsstrahlung heat- 
ing rates should influence laser plasma interaction ex- 
periments. For example, the coUisional heating of long 
preformed plasmas by high intensity short laser pulses 
accompanies several important applications, including X- 
ray lasers, laser wake-field particle accelerators, and the 
high-harmonics generation. Moderate ion densities in the 
plasma waveguide (n^ = 10^^ — lO'^^ cm~'^) should make 
it easier to satisfy condition Q. High intensity picosec- 
ond laser pulses would allow investigation of the heating 
rate dependence on ve/vtb over the wide range of val- 
ues. Atomic processes involved in X-ray lasers are very 
sensitive to plasma temperature and could provide a di- 
agnostic method for the inverse bremsstrahlung heating 
rates in addition to other diagnostics techniques such as 
the Thomson scattering and absorption measurements. 
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